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I. INTRODUCTION 


Among methods allowing for hyperthermia in biological tissues, that based on the radiofrequency (RF) 
absorption by gold nanoparticles (GNPs) got in the last years much attention, due to its important advantage 
- a deep penetration of the RF radiation into biological tissue (~30 cm at frequencies ~10 MHz). Details 
and useful references can be found in the review by Corr et al. (2012bI. Moran et al . I (2009) in their 


experiments carried out at 13.56 MHz revealed that the heating rate (HR) depends on the particle size, 
namely that smaller particles are heated faster, presumably because of their higher resistivity. However, 
the exact physical basis of heat generation by nanoparticles remained so far unclarified. This is a serious 
obstacle for grasping the HR / particle size relation and hence for making “intelligent choice” of the GNP 
size, tailored to the particular use. For example, one may wish to assure an enhanced specific absorption 
rate (SAR) to heat the GNP to necessary temperatures at minimum released power, avoiding at the same 
time an unnecessary heating of normal tissues during hyperthermia. 

According to standard definition, SAR = C(AT/At), where C is the specific heat and AT/At is the 
HR (temperature / time). The size dependence of the HR for the GNPs used by Moran et al. ( |2009 ) is 
presented in Fig. [l] As is seen, the HR attains maximal values, among the nanoparticles studied, in the 


smallest ones, of the 5 nm size. Along with Moran et al. (20091, also the observations of Kruse et al. (2011); 
Raoof and Curley] ( |2011 ) imply that the smaller the size of GNPs (for a given volume fraction of gold), the 
higher their HR: among experimentally studied GNPs with diverse sizes, those of ~5 nm diameter exhibit 


the highest one. Corr et al. (2012a) strengthened this conclusion, having stated that the GNPs only of size 


10 nm or less could be heated by the electromagnetic RF field (again, the highest experimentally observed 
HR corresponded to the GNPs of 5 nm diameter). 

In this paper, a physical model of the size effect in heat generation and reduction of heat generation 
in aggregated GNPs is proposed, in which the longitudinal acoustic vibrational modes (LAYMs) play an 


important role. None of the theoretical models reviewed critically by Corr et al. (2012b) took these latter 
into account, whereas we will show that the inclusion of acoustic modes enables to explain a number of 
experimental results. Besides, our findings make it possible to obtain an optimal size of GNPs providing 
maximum SAR. 

The manuscript is organised as follows. Sec. [TT] shows that in large GNPs (of the diameter exceeding 
162 nm) the heating by RF radiation is possible owing to uncertainty in the Fermi electron’s momentum; 
however, the use of this mechanism in hyperthermia cannot be but very limited. For practical treatment, 
much smaller GNPs need to be - and indeed are - used, revealing a pronounced size dependence of the 
heating efficiency. Possible physical foundations of the related heating mechanism, along with explanation 
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FIG. 1 Size-dependent heating of 
gold nanoparticles (diameters indi¬ 
cated) in the RF held of the frequency 
13.56 MHz. Reproduced with permis¬ 
sion from Moran et al. (20091. 































3 


of the observed size effect, are discussed in Sec. Ill In Sec. IIVI on the basis of conclusions drawn from 


Sec. |III[ the experimentally observed reduced heat production in aggregated GNPs is explained. Sec. [V] 
addresses possible ways of further increase of the HP in GNPs. The general discussion of results is offered 
in Sec. m and the conclusion drawn in Sec. m 


II. SETTING THE STAGE; CONDITIONS FOR THE RF ABSORPTION 

First we discuss an issue of whether the phonons are at all needed for the RF absorption, and what 
are the conditions for the latter to happen in a “conventional” way. The key feature in absorption of 
electromagnetic waves by nanoparticles is quantization of energy levels, which become denser as the particle 
size increases. Satisfying the conditions of conservation of both energy and momentum on absorption is not, 
a priori , possible without additional considerations, because the curvature/slope of dispersion relations are 
very different for electrons and photons. We consider first the case when no other players enter the picture, 
and assume the matching conditions to be satisfied within the uncertainty relation. To be specific, we take 
the GNP to be spherical, of diameter D\ the RF of 13.56 MHz, like in the above experiment, corresponds 
to photon energy hv = 5.6-10 -8 eV and momentum p p h ~ 3-10 -30 g-cm-s -1 . 

The skin depth in gold at 13.56 MHz is about 20 /rrn, i.e., by orders of magnitude exceeding the GNP 
size. One can therefore assume that the electric field penetrates the volume and is of the same strength 
throughout the particle. 

In metal particles as small as ~5 nm, and even for much larger ones, the quantization (discreteness) of 
the electron energy spectrum is a crucial factor shaping their properties. Speaking specifically about gold, 
one can advance far enough on the basis of the free-electron model, assuming one electron per atom and 
estimating the corresponding Fermi energy. 

Assume that the absorption happens due to transitions between discrete levels induced by quantum 
confinement; the free-electron dispersion relation E = p 2 /2m implies that the energy step A E = E p h must 
be accompanied by the correction of momentum, A pp ss Ep^m/pp = hv/v f (see Fig. [2]). For the relevant 
values of hv = 5.6T0 -8 eV and the Fermi velocity of gold as in the free-electron model (up ~ 1.4-10 6 m-s -1 ), 
the momentum mismatch amounts to A pp ~ 6.4-1CG 28 g-cm-s _1 , two orders of magnitude larger than what 
the photon momentum (see above) could bring about. However, the small size of GNP may result in 
appreciable Heisenberg’s uncertainty of electron momenta in them and thus “smear out” the discussed 
mismatch condition. Let us make some estimates to this end. 


Following the pioneering works by (Kubo 1962 19771, we remind a relationship between the diameter 
D and the level spacing A E. Assume that the electron energy level E' is the first vacant one above the 
Fermi level Ep. The level spacing AE e i = E' — Ep and the number of atoms N a in the GNP are related as 
follows (Kubo 1962 the derivation of the Kubo’s formula is given in Appendix |A|) : 


AE el » (4/3 )-Ep/N a , or N a « (4/3)-Ep/AE el ; 


D : 


8 m a Ep 


7T p A E e 


1/3 


(1) 



FIG. 2 Scheme of absorption of RF photon 
by the Fermi electron without involvement of 
LAVM. Dashed curve indicates a transition of 
the Fermi electron excited by the RF photon 
with energy hv. 
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FIG. 3 Schemes of “direct” absorption of RF photons hv by electrons at the Fermi level Ep , for different 
energy discretisation steps A E 1 depending on the particle size, and for different step multiplicities a (see 
text and Table [I] for details). The uncertainty of the electron’s momentum is hinted by elliptical dots of 
different width. 


where m a is the atomic mass of gold, m a = 197 in atomic units, or « 3.27-10 -22 g; p = 19.3 g-crn -3 
is the density of gold, Ep = 5.52 eV. For AE e \ = hv = 5.6T0 -5 meV, D = 162.1 nm. Remind that 
these are estimates for the smallest GNP size which allows the RF absorption at 13.56 MHz, from the 
energy level spacing arguments. Now, we turn to the momentum uncertainty arguments and assume 
for simplicity that the direction along the nanoparticle’s diameter is collinear with the momentum A pp 
and the Fermi electron’s momentum pp. Then the uncertainty in the Fermi electron’s momentum in this 
direction equals A p ~ h/(2it-D). At D = 162.1 nm, A p ~ 6.5T0 -23 g-cm-s _1 , that largely exceeds 
App ss 6.4-10~ 28 g-cm-s _1 . Consequently, the momentum conservation law can be helped by such uncer¬ 
tainty, and the RF photon absorbed. We keep in mind, however, that the minimal GNP size for the direct 
RF absorption must be at least «162.1 nm. Obviously as the particle size grows, the distance between 
quantified levels shrinks, so that an absorption of a given photon hv may bring about an excitation across 
a (an integer number) interlevel distances. This situation is schematically shown in Fig. [3j which implies 
a parabolic dispersion law E = p 2 /2m for different confinements (particle sizes), in arbitrary scale. The 
effect of the momentum uncertainty is indicated by dots of different widths placed at the initial and fi¬ 
nal levels. Some numerical estimates for the energy spacings and momentum uncertainties, relevant for 
the absorption of quanta at v = 13.56 MHz (hv = 5.6-10 -5 meV) are summarized in Table [I] An issue 
of matrix elements for such transitions, some of which might be suppressed in the dipole approximation 
due to symmetry considerations, is left aside. The numerical relation between the energy split multiplicity 
a, given by AE e i = hv/a , and the diameter of particle (assuming the latter spherical) from Eq. ([I]) is 
as follows: 


D^a 1 ' 3 - [(8/tt). (m a /p)-(Ep/hv)] 1 ^ 3 = a 1/3 -162 .1 [nm]. 


( 2 ) 


TABLE I Parameters of GNPs with level spacings hv/a (a=l, 2, 4, 10 taken as examples) able to absorb 
the 13.56 MHz photons thanks to the uncertainties in the Fermi electrons’ momenta 


Number of 
level spacings a 

Level spacing 
hv/a (10 -5 meV) 

Diameter 

D a (nm) 

Uncertainty in Fermi 
electron’s momentum A p a 
(10~ 23 g-cnns -1 ) 

1 

5.6 

162.1 

6.5 

2 

2.8 

204.2 

5.2 

4 

1.4 

257.3 

4.1 

10 

0.6 

349.2 

3.0 
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The last column of Table [I] specifies the corresponding momentum uncertainty, Ii/(2tt-D). In all the cases 
considered (until entering much larger particle sizes than those covered by Table |l]) , the uncertainty in the 
Fermi electron’s momentum A p a exceeds by far the “adjustment” A pp k, 6.4-10 -28 g-cm-s _1 , discussed 
above as a reference value to absorb the momentum mismatch required to promote a Fermi electron onto 
the energy level (Ep + hv). Thus, the momentum conservation law can be fulfilled, and the 13.56 MHz 
photons absorbed, for GNPs larger than ~162 nm. Consequently the latter, as it seems, can apparently 
be RF-heated “directly”, i.e., without involvement of LAVMs, just by exciting the Fermi electrons to the 
available energy levels above the Fermi level. On the contrary, in GNPs smaller than ~162 nm the separation 
between adjacent levels exceeds the energy of the RF quantum, and at “working” particle sizes which in 
fact show an enhanced absorption (~5 nm) the level splitting becomes prohibitively large. An involvement 
of many photons in a single electron excitation is likely to be a very rare process. Therefore, we have to 
bring another physical mechanism into consideration. Longitudinal acoustic phonons seem to be able to 
intervene with energies of the “correct” order of magnitude. The prerequisites of their practical involvement 
are discussed in the following. 


III. LAVM-ASSISTED ABSORPTION OF RF PHOTON BY A FERMI ELECTRON 


The observed size dependence of the HR can be explained by a mechanism of the GNPs’ heating that 
attributes a crucial role in the absorption of a RF photon to LAVM. To make reliable estimates, we need 


to know something about elastic properties of nanoparticles. A theoretical work by D.Y. Sun et al. (2001) 


demonstrated that the metallic nanoparticles retain the bulklike core region. Experimentally, J. Sun et oT] 
(2014) have recently found that Ag nanoparticles can be deformed like liquid droplets but remain highly 


crystalline in the interior. 

As was already admitted above, we assume the GNP to be spherical; we remind that the RF of 13.56 MHz 
corresponds to photon energy kv = 5.6T0~ 8 eV and momentum p p h~34CV 30 g-cm-s _1 . The absorption 
happens due to transitions between discrete levels induced by quantum confinement, as was elaborated above 
and schematically shown in Fig. [3j 

From now on, we turn to discussing a case of nanoparticles too small for a RF phonon energy to 
“bridge the gap” between largely split quantified levels. We bring into consideration a scheme of absorption 
of a RF photon whereby a Fermi-level electron excitation is helped by an involvement of a LAVM. A 
“classical” view onto the electron-phonon interaction (otherwise straightforwardly grasped in terms of energy 
and momentum exchange between the corresponding quasiparticles) is that the electrons are driven by, or 
themselves contribute to, the fluctuating electric field due to compression / dilation of the electron density 
in the course of lattice vibrations. Obviously, only the longitudinal vibration mode can be “useful” in this 
sense. Its wavevector moreover must be reasonably far from the Brillouin zone (BZ) center, where the 
dispersion starts to bend “flat” and to yield high density of modesQ 

A Fermi electron may absorb the energies of both the RF photon and the LAVM. Fig. [4] shows how an 
inclusion of momentum n vm -q and energy n vm AE vm of such LAVM ( n vm numbers the energy steps in the 
vibration modes) into the combined absorption makes possible to satisfy both conservation laws. In the 
absorption event, the GNP borrows energy from its LAVMs system, adds it to that of the RF quantum and 
excites an electron to the vacant level beyond the Ep: 


Ep + A E vm + hv = (p F + Q + P P h) 2 /2m , 


( 3 ) 


where Ep/pp , AE vm /q and hv/p v p are the energy / momentum of the free electron, the LAVM and photon, 
correspondingly (see Fig. [4]}, and m is the electron’s mass. We note, for the sake of a later reference, that 
an electron may well be excited across m e i > 1 steps of its discrete spectrum, i.e., onto Ep + m e \AE e \, 
that is however not explicitly depicted in Fig. [4] Once excited, the propagating mobile electron very likely 
will be trapped (relaxed) before reaching the GNP surface (see Subsec. III.C on the mean free path issues), 
releasing more energy to the LAVMs pool than was borrowed beforehand. The net effect of that is the 
heating of the GNP. 


X A numerical estimate can be drawn from (Singh and Ali 2013), who studied the phonon dispersions in amorphous 
metals. For gold, the first maxiumum of the io(q) dispersion occurs at g~1.5 A -1 whereas up to q~1.0 A -1 the 


dispersion remains reasonably linear. 
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FIG. 4 Scheme of absorption of RF photon 
by the Fermi electron assisted by an absorp¬ 
tion of a LAVM (in arbitrary scale). The dis¬ 
persion laws for photons Aph^Ppiu electrons 
E(p)~p 2 /2m and LAVMs E vm ~(n vm -q) are 
shown. See text for details. 


In Fig. |4j a point on the Fermi sphere (identified by momentum pp and energy Ep) serves as an origin of 
the LAVM dispersion branch, on which the momentum n vm ■ q and the mode energy A E vm stand for some 
“representative” mode. Let us give some numerical estimates to these parameters. 

The free-electron model for bulk gold yields the Fermi radius kp = 1.20T0 10 m _1 , quite close to the 
mean experimental estimates for its slightly non-spherical belly (Coleridge and Templeton 1982). The 


corresponding free-electron Fermi energy is Ep = 5.52 eV; the Kubo formula for the level splitting at the 
Fermi energy due to spatial confinement A E vm yields, for the diameter D = 5 nm, A E vm = 1.91 meV, 
and the corresponding increment of the electron’s momentum to bring it to the first vacant level is A kp = 
2.08-10 6 nr -1 « 0.017 of the kp. As RF hv = 5.6T0 -5 meV<C A E vm , the larger part of the latter gap 
has to be overcome, according to out hypothesis, by borrowing energy from LAVM. From the velocity of 
sound in gold (3240 m-s _1 longitudinal, 1200 m-s -1 transversal), the highest phonon energy extrapolated 
onto the nearest BZ boundary, i.e., the L point, is ~28.5 meV (Lp) and ^10.5 meV (Lt)- Corresponding 
experimental frequencies at L are bended downwards as expected, to, correspondingly, 19.1 and 7.7 meV 
(as cited by Dal Corso 2013 see also his first-principles calculation of vibration spectrum). Consequently, 
the energy matching to A E vm can be realized via an interaction with an acoustic phonon whose momentum 
is within ~7 — 18% of the BZ radius. Inversely, an interaction with a phonon close to the BZ periphery 
may, in principle, promote an electron onto a higher (up to the 3th or 4tli) energy state beyond the Fermi 
level. Two observations, however, need to be stressed in this relation. First, the numerical relations in 
gold (see above) are such that matching the energies of electron and phonon is impossible if their momenta 
stay collinear, as Fig. [ 4 ] implies. A more realistic scenario, discussed below, comprises kp staying at some 
angle to the phonon momentum q. The second observation imposes the quantisation of phonon q values in 
a nanoparticle as q = h/L , where L is a length of the chain of ions along which the LAVM propagates. As 
an estimate of the order of magnitude, the relevant length is between the nanopaticle’s diameter (5 nm) and 
circumference, hence the q step is A q = [1- • -2tt\/L = [0.2- • -1.26JT0 9 m -1 , i.e., ~[2-ll]% of kp. 

In the following, we’ll often refer to circular contour / path over which the compression mode propagates, 
as this is a simple yet realistic model case. Although the natural diversity of nanoparticle sizes and shapes 
makes a faithful simulation difficult, the circular path has a virtue of being the longest one in a “round” 
particle, hence hosting the maximal number of modes, densely distributed in the q space. Consequently, the 
quantified energies of modes on a circular contour make a denser spectrum than those on any other path; 
within a given energy interval, more individual modes can be found and used for borrowing energy to an 
electron. It will be argued below that longitudinal acoustic phonon cannot propagate strictly on the surface, 
but rather at a (small) depth. 

We discuss now the impact of the electron and phonon momenta being non-collinear, but, in any case, 
coplanar. A more realistic scheme in this sense than that of Fig. [4] is shown in Fig. [5] where however the 
energy and the momentum of the RF photon are neglected, by force of relations hv <C E vm ; p p p <C n vm Q', 
Pph pp. Fig. [5] depicts therefore a possibility of absorption, by electrons at the Fermi surface, of the 
LAVMs with energies A E vm , 2A E vm , 3A E vm , and corresponding momenta, within the same nanoparticle 
of a given size. It is (arbitrarily, just for the sake of simplifying the figure) implied that m e \ = n vm and 
hence AF e ] = A E vmi therefore the augmentation of the electron energy on absorption happens in portions 
of AE vm . 
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FIG. 5 Same as in Fig. [1J however assuming 
m e i = n vm and A.E e \ = A E vm , for the gen¬ 
eral case of the electron momentum pp and the 
phonon momentum q being coplanar. For sim¬ 
plicity, the energy and the momentum of the 
RF photon are neglected. See text for details. 


A more detailed projection of the momenta matching, involving also the RF photon momentum, is shown 
in Fig. [6] Possible excitations from an initial state pp of energy Ep via an absorption of a phonon q and a 
RF photon p p h end up in a state with the energy E = Ep + A E vm + hv. The allowed “chained” pp- ■ -q- ■ -p p h 
vectors fall within a body of revolution around the fixed p p h||i?x direction, limited on the left by the cone 
of [side pp / aperture 2/3] and on the right - by the spherical (radius \pp + q + p p h|) cape built on top of 
the cylinder of the height p p p. The conical and cylindrical surface parts are connected by the intermediate 
conical belt of the width q. From Fig. |6][a), the maximal angle /3 the pp vector may build to q is 


/? = arccos 


(V2 m e E - q) 2 - 

2(Pph\PF) 


(4) 


On relaxation of the excited electron back to Ep , the extra energy E' = AE vm + hv is released into the 
vibration pool, exciting a phonon with \<f \ > |cf|. hv is therefore the net gain in energy. 




photon momentum p p ^, neglected in Fig. [5j is now retained. Circles centered at the origin of the coordinate 
system, of radii \pp \ and \pp + q + p p h |, are cross-sections of the electron energy paraboloid at E = Ep and 
E = Ep + A E vm + hv. (a): Some examples of matching momenta and energies in the process of electron 
excitation from pp due to absorption of phonon q and RF photon p p p. ( b ): For one of pp vectors from 
panel (a), excitation channels are followed by possible relaxation channels into the initial state via release 
of phonon if. Note that \cf\ > \q\. 
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FIG. 7 VDOS in AU 959 nanocrystal, re¬ 
produced with our modification with per¬ 
mission from (Sun et al. 2001). Curve 1: 
contribution from the surface shell; curve 
2: from the transition shell, curve 3: from 
the core region; curve 4: our approximate 
extraction of the VDOS for LAVMs in the 
core region, yielding a peak centered at 
3.7 THz (15.3 meV), with full width at half 
maximum (FWHM) ~1 THz (~4.14 meV). 


In the following, we consider the LAVM within a GNP as a compression wave propagating along a closed 
chain of atoms, in the spirit of the Born - von Karman cyclic boundary conditions. For this analysis, the 


work by D.Y. (Sun et al. 2001) is useful which deals with elastic properties and vibrational density of 


states (VDOS) in GNPs, indicating notably three structure elements of a nanoparticle: the surface shell (of 
« 2 A thickness), the transition shell (of ~ 3 A thickness beneath) and the core region. Fig. [t] depicts the 
corresponding decomposition of VDOS for the case of 959-atom GNP, of the diameter ss 3.2 nm. In the core 
and transition shell, the distribution of VDOS clearly shows the peaks at « 3.7 THz (E vm ss 15.3 meV), 
resembling the peak in the VDOS of bulk gold in the experimental work by (M unoz et a7| [2013|, which we 
attribute as being due to LAVMs (see Fig. [ 8 ]). The VDOS of the surface shell in Fig. [T] does also indicate a 
feature at the corresponding energy, that we extracted from the background as the curve 4 in Fig. [7| 

The VDOS maximum under discussion is strongest in the core region, which is therefore expected to 
contribute mostly to the RF absorption. According to D.Y. (Sun et al. 2001), the core region is bulklike, 


therefore while discussing its contribution to the RF absorption, one can rely on the density and sound ve¬ 
locity values for crystalline gold. This is consistent with experimental observations on the gold’s homologue, 


silver, by J. (Sun et al. 2014) who have found that its nanoparticles remain highly crystalline in the interior. 


In the following discussion, we assume that the “useful” LAVMs propagate at the depth 5 = 0.7 nm, i.e., 
along the closed contour entirely within the bulklike core region. 



FIG. 8 Experimental phonon DOS curves 
for bulk pure Au and Au0.97Fe0.03 from 
the inelastic neutron scattering measure¬ 
ments; reproduced with permission from 
( Munoz et al . | 20131. Dashed line is our 
approximate extraction of the VDOS curve 
for LAVMs (the FWHM range is from 
14.6 to 18.4 meV; the peak position is at 
16.8 meV). 
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A. Energy conservation 


We discuss now some quantization relations for electrons and LAVMs following from the GNP geometry. 
The energy conservation condition for the absorption of a RF photon with the energy hv is as follows: 


m eF AE e \ — n vm -AE vm T hv 


Tl vm 


•A E vm . 


(5) 


Here, m e \ is the number of steps (gaps) between the quantized electron levels, and n vm the number of 
vibration quanta helping a RF photon to get absorbed. According to the Kubo’s formula (Kubo [1962 


1977) - see also Eq. (A4) in the Appendix [Aj - the step in the electron energy levels A2? e i depends on the 
number of gold atoms N a . The condition (Eq. [5]) that m e \ energy steps must embrace n vm phonon energies 


takes the form 


4 E F ^ n vm v L h 
3 TOel Na ~ 7 t(D - 2(5) ’ 


( 6 ) 


where L = n(D — 26) is the length of the closed contour at the depth 6 under the GNP surface, and the 
(longitudinal) sound velocity v F relates frequency to wave vector. Further on, assuming that the density of 
gold in the surface shell and the transition shell is close to that in the core region, i.e. in bulk gold, the 
number of atoms in GNP can be expressed via density of gold p, atomic mass m a and the volume V (or 
diameter D) of the particle: 


* jl d ». 

6 m a 


m„ 


Taken together with Eq. ([6]), this yields the depressed cubic equation on D: 

+ p D + q = 0 , 


(7) 

( 8 ) 


with p 


8 E F m a ( TOel 
V F hp V n vm 


q = 


\QE f m a 6 
vl hp 



(9) 


For p < 0 and q > 0, Eq. ([8]) always has one negative root and either two (possibly degenerate) or none 
positive ones. The two positive roots of Eq. ([8]) can be expressed as follows: 


D k = —2 




arccos 


q 

2 




( 10 ) 


where k = 1,2. Note that the parameters p and q contain, along with constants depending on the properties 
of gold, the trial numbers 6 , ?n e i and n vm ; recall that <5 is the depth of propagation of LAVM, and n vm is the 
number of vibration quanta matching within the electronic excitation. The solutions D k can be expressed 
in terms of 6 , m e \ and n vm as follows (D k and 6 measured in nm): 


D k = 



+ 


2n k' 
~ 3 ~. ' 


( 11 ) 


Of two positive roots given by Eq. we retain the practically relevant largest value of D k , corresponding 
to k = 1. Note that certain combinations (m e i, n vm ) yield the argument of arccosine > 1 and hence no 
solution. The roots are densely yet unevenly distributed, as is shown by Fig. [9] Note that the total number 
of solutions increases, as expected, with the number of trial (?n e i, n vm ) combinations, however the profile of 
the root density, with its narrow maximum around « 7 nm, remains stable up to the upper ~ 25% of the 
total span of D values. For this reason, and moreover since a non-ideal equidistanteness of both electrons’ 
and LAVMs’ levels would eventually “detune” the criteria set by Eqs. for large (TO e i, n V m)> we set, 

from now on, an arbitrary limit (= 10) on the maximal values of the latter. 

We’ll see below that this remarkable “selectivity” of GNP sizes D with respect to their capacity to satisfy 
the energy conservation relations will eventually manifest itself in the distribution of the HR. The latter 
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FIG. 9 The density of discrete solutions Dk =1 
after Eq. ( [II] ), smeared into a continuous dis¬ 
tribution with the half-width parameter of 
1 nm, for better visibility. <5 is fixed to 0.7 nm. 
The cases (m e i, n vm )<10 (red dots, refer to the 
left scale) and (m e i,n vm )<20 (blue continuous 
line, refer to the right scale) are compared. See 
text for details. 


ought to be influenced, however, by the next element entering our discussion, namely, the availability of 
“useful” LAVMs that can be induced in the particles of selected sizes. “Useful” means the modes whose 
energies are multiples of n vm -AE vm = n vm -vi J h/['K( y D — 25)], n vm being selected by the commensurability 
of the vibration energies with electron excitation ones, in the spirit of Eq. (§. Obviously, for the optimal 
heating it is essential to find many such vibrational modes within the FWHM of the LAVM-related VDOS 
peak of gold (cf. Figs. [7|8] and the related discussion), the number we’ll refer to as -/Vfwhm in the following. 


B. Momentum conservation 

In addition to the energy conservation equation ([5]), one should take into account the momentum con¬ 
servation condition. We’ll specify it for the case of LAVM propagating along the circular contour of the 
diameter D — 25, to which the phonon momentum n vm q is tangential - see Fig. |10| Anywhere on the contour, 
a Fermi electron with the momentum pp can intervene to bring about an absorption of a RF photon (we 
neglect the effect of the momentum and energy of the latter onto the resulting conservation relation, as was 
already argued before). 

The situation arbitrarily shown in Fig. [lO] assumes the momentum of the Fermi electron to point inwards 
the GNP; on absorbing a phonon, such electron would transverse the particle along the path l and release 
the phonon “on the other side” of the contour. A different possibility would be the pp pointing outwards; on 
the nanoparticle’s surface such electron would be either elastically reflected inwards, or emit a quantum and 
be “lost” for the GNP heating mechanism we discuss. The probability relation for such elastic / inelastic 


FIG. 10 Scheme of RF photon absorption on a cir¬ 
cular contour inside a nanoparticle; see text for dis- 


S 



cussion. 













11 


events at the surface is a priori difficult to estimate. Anyway, the elastically scattered electron will almost 
“instantly” regain the contour and follow the path of the inwards-moving electron, as depicted in Fig. [Toj 
To specify what “instantly” means, we note that the perturbation of the potential inside the particle, “felt” 
by an electron, propagates, at most, with the longitudinal sound velocity in gold, i.e., v l = 3.23-10 5 cm-s _1 , 
whereas the Fermi velocity in gold is three orders of magnitude larger: the free-electron model with the 
electron density corresponding to that of monovalent gold yields vp = 1.4-10 8 cnrs ' 1 (cf. Sec. [n| see also 


Ashcroft and Mermin 19761). The “flight time” for the electron to reach the surface of the nanoparticle and 


get back to the contour is therefore ~ 10~ i5 s. 

A simple geometric argument illustrated by Fig. 10 (for the ideal case of planar circular contour) shows 
that the momentum of the absorbed phonon n vm q can “ride” the electron across the particle and be released 
on the other side of the contour, that we’ll refer to as relaxation. In fact, simultaneously released are 
the Fermi momentum pp and the phonon momentum n vm (f, both at angles with their respective “pre¬ 
absorption” values but exactly preserving the corresponding moduli, under the condition that the phonon 
is emitted along the contour at the electron’s contact with the latter on arrival. Specifically in Fig. [lOj the 
transferred momentum s, a sum of the electron pp and phonon n vm q momenta are related to the electron 
emission angle 7 as follows: 


Pp= s 2 + ( n vm -q ) 2 - 2s(n vm -q) cos 7 , 

i.e., the angle the momentum of the excited electron makes to the contour is 


7 = arccos 


s 2 + ( n vm -q ) 2 -p| 
2 s(n vm -q) 


( 12 ) 


(13) 


Fig.[lO]implies moreover that the momentum of a LAVM phonon, via the interaction with an electron, may be 
“reinforced” by that of RF photon. The “parallel” (along the contour) component of the latter is then added 
to n vm -q and transferred (by mediation of the excited electron) to the relaxation point. The momentum of 
then released phonon is n vm '<f , with |<f | > \q\, even if the net increase of the phonon momentum is quite 
small, with respect to electron and phonon counterparts: hv <C A E vm ; p p h <C nq , p p h -C Pf- 

The magnitude of s in Eq. (13) is s = y/2m(Ep + n vla AE vm ). On expressing the magnitude of the 
LAVM momentum q and the energy step via the nanoparticle size D and the velocity of (longitudinal) 
sound ul, 


q = h/[ir(D - 2d)]; A E vm = vph/[n(D - 2d)], 
we arrive at an expression for 7 in terms of n vm and (D — 2 d): 

?l V m h 


7 = arccos ■ 


mv l + 


2t t{D - 2d) 


2 2 m n vm vl h 

Pf + 


7 t(D - 2d) 


(14) 


(15) 


An excited electron would traverse the nanoparticle (along the chord l in Fig. 10 1 and release energy on 
entering the cyclic contour again, whereby a vibration mode with momentum n vin (f would be induced. 
Technically this might happen as a consequence of an electric field being suddenly created at the “arrival 
point” of the electron on the contour; the Coulombic attraction of the ions would trigger the compression 
wave to run along the contour. The electron path (chord) is related to the nanoparticle parameters n vm , D , 
d as follows: 


l = (D-2S) cos- 7 ) 
= (D — 2d) sin arccos < 


mv l 


2tt(D — 2d) 


2 2m n vm vph 
Pf ' 


n(D - 2d) _ 


(16) 


For relevant values of n vm and (D— 2d), the argument of arccosine varies within ~ 0.006—0.187. Consequently 
the sine in the above formula stays within 0.982 — 0.999, i.e., I ~ (D — 2d) with, at most, ~ 2% accuracy, 
therefore the excited electron must transverse the nanoparticle roughly along the latter’s axis. On neglecting 
the discreteness on the circle of LAVM propagation (of the D — 26 diameter), each its point may serve as 
the “source” of the excited electron. The regions of absorption of RF photons in the nanoparticle are the 
rings of (D — 2 d) diameter. 
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C. Electron free path 


The probability for an electron to transverse the particle, i.e., to escape being scattered along the path of 
the length l and to reach the “opposite” point on the cyclic contour, equals exp(—Z/Zo)~ exp[— (D — 25) /Iq], 
where Iq is the mean free path of an electron in the nanoparticle of the size D (see below). 

It seems logical that the HR of a nanoparticle of size D is proportional to the following factors: (i) 
the energy of absorbed RF photons, (ii) the length of the contour ^r(D — 25) (as the absorption may 
occur in any point thereupon); (in) the summary number E,; ^fwhm of matching opportunities for the 
multiple energy step n vm -AE vm to fall within the energy range of FWHM vibration modes [the estimates 

or 12.86 — 16.54 meV according to D.Y. (Sun 


are 14 .60 — 18.40 meV according to (|Munoz et al. 
et al. 


2013) 


2001 )]|^] and (iv) the probability for an excited electron to undergo relaxation on the contour, i.e. 


exp[— (D - 2 S)/l 0 \: 


HR - 


hv-2 7T (D - 26) Y n fwhm exp[-(T> - 26)/l 0 \ 


(17) 


Here, the factor 2 accounts for two senses of propagation of vibration modes along the closed contour, and 
the index *, without further elaborating for the moment, identifies absorption / relaxation “events” likely to 
contribute to the HR. This summation should, at least, take into account different combinations (m e i, n vm ), 
within the global limitation <10 imposed on these parameters, which retain m e i/n vm and thus correspond 
to the same solution D of Eq. but whose values IVfwhm are different. In the next subsection, we’ll 

adress an issue of multiple contours which may participate in the absorption and the relaxation events, in 
which relation the summation will be further explained. For the time being, we retain the general structure 
of the expected expression for the HR, and try to specify the relevant value of Iq ■ The upper limit on it is the 
mean free path in bulk gold, ~ 35 nm; in nanoparticles, due to irregularities of internal structure, Iq ought 
to be much shorter. A priori, the smaller the size, the more likely the crystal lattice is distorted, hence the 
smaller the mean free path. Estimates for some “working” value for small enough GNP sizes, which were 
covered by the study by (Moran et al. 20091, can be gained from their experimental data reproduced in 


Fig. 0 with the help of insight given by Eq. (|17|. 

Fig. [l] reveals a tendency towards the “saturation” of curves with an increase of the gold volume fraction. 
This reflects the loss of efficiency of heating the nanoparticles following their aggregation: the higher the 
GNP concentration, the higher the probability of their agglutination. Therefore the “net” Iq values are more 
safely to extract from the slope of curves d(HR)/d(volume fraction) near the origin. Of interest for us is 
the value of Iq for GNPs with D 1 = 5 nm and D 2 = 10 nm; the corresponding slopes are almost identical. 


Eq. (17) yields the following relation for the curves corresponding to these values of D: 


HR 2 ~ D 2 -26 Ef 2) <hm 


HR, 


D , - 25 


E 


(-Di) 7\r0) 

iV FWHM 


exp 


D 1 — D$ 


(18) 


As mentioned above, each sum implies all relevant absorption / relaxation events within the particle of the 
corresponding size. Practical calculations show that, in order to estimate the HR 2 /HR 1 ratio from the data 
of Fig. [I] with the accuracy of « 10%, it suffices to retain in Eq. (18) the leading term of each sum. Then, 


for the range of diameters ~ 5 — 10 nm, l 0 can be expressed as follows: 


k = (£>1 - D 2 ) 


In 


HR. 


2 Di 25 A fp^ HM \ 


HRi D 2 - 25 at (2) / 


(19) 


The effective identity of the d(HR)/cZ(volume fraction) slopes characterizing the GNPs of D\ = 5 nm and 
D 2 = 10 nm implies HR 2 /HR 1 WI in Eq. (18). For the contour 5 = 0.7 nm deep, (D\~25) / (D 2 — 25) ~ 0.42. 


The relation -Npwhm/-^fwhm follows from the straightforward counting of LAVMs that make a discrete 
spectrum on a circuar contour of the (D — 26) size: how many of modes will fall within the FWHM of the 
longitudinal acoustic peak of gold. In principle, -/Vfwhm steadily increases with size, but the subtlety is that 


r(2) 


2 The sum JA IVp^ HM takes into consideration contributions from the absorption of RF photons at phononic momenta 
q, 2q, iq etc., in the spirit of Fig. [5j for fixed m e i/n vm and <5. 
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not all steps in vibration energy are compatible with the energy conservation criteria Eq. 00 Appendix 
b] explains this situation and argues that -Np^ HM = 1, « 7.5, hence -Npwhm/-^fwhm « 0.13, and 

—5 nm/ln(0.055) = 1.72 nm. 

.. as 


from Eq. (191 Iq 


Coming back to the discussion at this subsection’s opening, we can anticipate that the final HR, 
function of GNP size, will be the interplay (multiplication) of three tendencies. The first one is the availability 
of “good” D values which can contribute at all; they group around ~7 nm and rapidly become scarse at 
larger sizes. The second tendency is the number of LAVMs within the “good” energy interval (given by the 
elastic properties of gold); this number essentially grows with D. The third effect is an exponential cutting 
of the HR at the characteristic length much shorter than the electron mean free path in the bulk gold. 
In total, the last tendency shifts the maximum of HR(D) a bit more to the left from the abovementioned 
“primary” « 7 nm value than the second tendency shifts it to the right. 

One can infer that, would the data for D < 5 nm be available in Fig. |T| one could expect the initial 
slopes of the corresponding curves (of HR vs gold volume fraction) to go steeper. Correspondingly, the Iq 
values for such sizes would likely be less than 1.72 nm, and, on the scale of the GNP sizes, the maximum of 
the HR would occur at slightly smaller values than so far reported. 

We’ll see in the following how the final counting of HRs proceeds, which also takes into account a somehow 
delicate issue of possible “diversification” of the absorption / relaxation events. 


D. Case of multiple contours 

At a risk of attributing too much precision to a simple enough model, we would like to emphasize 
a possibility for the energy of excited electron to be returned to a “different” phonon than that it was 
originally borrowed from. Once N a and m e i are fixed for a GNP, the Eq. ([6]) can be, in principle, satisfied 
for various values of S' > 6, each S' being a different “depth” of a circular contour, assuming for simplicity 
a symmetric placement of the latter within the particle. The phonon energy at each contour is quantified, 
so that 


(D - 2 S')/{D - 25) = < m /n vm , (20) 

for integer n' vm . Obviously, the case n vm = 1 does not permit any inner contour, whereas n vm > 1 allows 
n' vnl = 1 ,.. ., n vm — 1. This is depicted in Fig. |TT] for n vm = 1 , 2 , 3, and summarized in Table [TT| In this way, 
the discretized phonon momentum values (integer) xh/L and the corresponding (momentum) xul energies 
must find their “resonance” counterparts on relaxation. This is possible if the “secondary” contour length is 
commensurate with the “primary” one. The relation is not exact, since the energy absorbed by an electron 
is that of phonon(s) plus (much smaller) RF quantum, whereas the “relaxation” energy (released by the 
electron) is just the phonon(s)’ one. However, it is helpful for counting different contributions. An essential 
observation is that the density of “resonance” modes decreases as the contour length shrinks. Still, the 
contributions from “inner” contours is not negligible; counting them, with different n vm /n' v m integer taken 
into account, would modify somehow the contributions of different GNP sizes to the HR. 

For each “event”, the “free path” of an excited electron U prior to its relaxation can be easily calculated 
in analogy to how it is shown, for several selected (n vm , n' vm ) combinations, in the last column of Table [TTJ 


TABLE II Characteristics of internal contours and the electron path transversing them for several small 
values of n vm . See text for details and refer to Fig. [Tl] 



< m 

CM 

1 

S 

CM 

1 

l/(D - 2 S) 

2 

1 

1/2 
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3 

1 

1/3 
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3 

2 

2/3 

l/6 e ; 5/6 f 

4 

1 

1/4 

3/8; 5/8 

4 

2 

2/4 

2/8; 6/8 

4 

3 

3/4 

1/8; 7/8 


Fig. 11 middle panel: a AB, b AC. Fig. 11 right panel: C AC, d AD , e AB, / AE. 
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FIG. 11 Schematically drawn cyclic contours within nanoparticles for three different values of the n vm 
parameter. “A” marks the excitation site for a Fermi electron on the outer contour, following the absorption 
of a RF photon and the LAVM. Other points indicate possible sites of relaxation of the excited electron on 
traversing the particle. 


in terms of the contour size. The values are listed in pairs, indicating that, as the electron traverses the 
GNP approximately along the diameter (by force of earlier envoked arguments), its “exit point” may occur 
on a close or on a remote point on a given internal contour. 

Alternatively, a situation is imaginable that the resonance criteria of Eqs. allow a modification of 

electron energy step on exchange with a phonon, on respecting the condition 


vh 


vh 


m e i-n(D — 25) m' el -n(D — 25') 


hence 


D — 25' 


TO e l 


( 21 ) 


Since the last relaton is < 1 (the secondary contour at depth 5' is deeper than the primary one), m' el > m e i; 
moreover, m' eX /m e \ must be an integer. 

Possible contributions to the HR, in the spirit of Eq. are summarized in Table [TTT| and depicted in 
Fig. 12 The summation in the fourth column of Table m selects the cases ( m e \ unchanged, n' vm variable) 
while that in the fifth column - the cases (m' el variable, n vm unchanged). 

The condition m' el < 10 combined with that of m' el /m e i to be integer does considerably restrict the 
amount of inner contours; hence so few non-zero entries in the fifth column of Table m The largest contri¬ 
butions (for D=6.65 and 9.80 nm) come out because, with their small initial m e \ value, the largest number 
of contours (= 4) could be generated. The summation over possible contours and transitions yields a global 



FIG. 12 Contributions of possible con¬ 
tours in GNPs of different sizes to their 
HR. Vertical lines mark the “individ¬ 
ual” values of -^fw HM ex p(— h/lo) for 
each “valid” contour within a “resonant” 
GNP size, whereby thin black lines stand 
for contributions of the {m' el variable, 
n vm = const) type, and thick orange 
lines - for contributions of the ( m e \ = 
const, n' vm variable) type. Light gray 
bars make a histogram of the abovemen- 
tioned contributions, summed up within 
2 nm steps of GNP diameter D. 
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TABLE III Contributions to the HR from summations over primary and corresponding secondary contours 
characterized by different (m e 1 , n vm ) values in gold nanoparticles of diameters D. See text for details. 






Ei^FWHM ex P( h/lo) 

m e i 

^vm 

D (nm) 

Ei — En; n 

,.< m <!0 Ei^Em'p^lel <<1<10 

1 

4 

2.44 

0.55 

0 

1 

3 

3.26 

0.76 

0 

2 

5 

3.74 

0.64 

0 

1 

2 

4.36 

1.10 

0 

4 

7 

4.75 

0.45 

0 

3 

5 

4.90 

0.50 

0 

2 

3 

5.22 

0.82 

0 

5 

7 

5.45 

0.72 

0 

3 

4 

5.61 

1.41 

0 

5 

6 

5.98 

0.92 

0 

1 

1 

6.65 

0.24 

4.62 

8 

7 

7.18 

0.42 

0 

7 

6 

7.26 

1.34 

0 

5 

4 

7.55 

0.87 

0 

9 

7 

7.67 

0.18 

0 

4 

3 

7.83 

0.79 

0.45 

10 

7 

8.14 

0.31 

0 

3 

2 

8.36 

1.51 

1.10 

5 

3 

8.86 

0.61 

0.38 

7 

4 

9.11 

1.24 

0 

2 

1 

9.80 

0.06 

3.13 

9 

4 

10.44 

0.94 

0 

7 

3 

10.65 

0.50 

0 

5 

2 

11.06 

1.10 

0.52 

8 

3 

11.45 

0.42 

0 

3 

1 

12.19 

0.02 

1.48 

10 

3 

12.90 

0.36 

0 

7 

2 

13.24 

0.93 

0 

4 

1 

14.21 

0.01 

1.00 

5 

1 

15.98 

0 

0.79 

6 

1 

17.58 

0 

0 

7 

1 

19.06 

0 

0 

8 

1 

20.43 

0 

0 

9 

1 

21.71 

0 

0 

10 

1 

22.93 

0 

0 


scan of the property JA -^fwhm ex P(~h/lo) which is expected to be proportional to the HR, as function of 


the GNP diameter. The results are shown in Fig. 12 in two ways: as a contribution (weighted by e li / l °) of 


each relaxation event, possible in a GNP of given size, and as cumulate effect of multiple events, grouped 
in a histogram with a fixed step in D. The first representation indicates that the contributions start from 
D = 2.44 nnr and drop down to zero beyond D = 15.98 nnr. The histogram representation seems more 
helpful in inspecting the “importance” of different GNP sizes for the HR. It seems that the HR efficiency 
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peaks around D ss 6 nm and rapidly decreases for smaller and larger diameters. 

To conclude the discussion about multiple contours, we point out that the excitation energy for an electron 
can also be borrowed from LAVM propagating along one of internal contours. However, the probability of 
such effects rapidly decreases with descending onto ever shorter contours characterized by ever sparser 
distribution of quantized resonance modes. 


IV. ON THE REDUCED HEAT PRODUCTION IN AGGREGATED GNPS 

An aggregation of n GNPs yields a larger particle, with the number of atoms and electrons increased by 
the factor of n. However, this won’t be normally accompanied by a formation of joint subsurface contour at 
the depth <5, as the core regions of different GNPs remain isolated from each other by their surface shells. 
Therefore the “optimal” conditions of RF photon absorption, given by Eqs. © and © , would be violated, 
and the HR of aggregated GNPs reduced. 

After Fig. |T| the HR eventually saturates, for all particle sizes, as function of gold volume fraction. For 
small GNPs this saturation occurs faster, for the apparent reason that, for a given gold volume fraction, 
smaller particles mean their higher concentration, and hence higher tendency for aggregation. 


V. FURTHER SUGGESTIONS FOR ENHANCING HEATING RATES IN GNPS 


In Sec. US it has been shown that in GNPs, the HR is enhanced thanks to large number of participating 
LAVMs. As heating of GNPs involves both the LAVMs and the Fermi electrons, one way of enhancing the 
HR ever further would be to increase the number of electrons absorbing the RF photons. This is possible 
through enhancing the local electronic DOS (EDOS) at the Fermi energy of gold, doping the latter with 


transition metal impurities, e.g., Ta or Fc see Fig. [13] and also (Moldosanov and Postnikov 2013 Postnikov 


and Moldosanov 2012). The main peaks of spin-split local EDOS of these impurities, according to first- 


principles calculations]^] are pinned at the Fermi level of gold, adding substantially to the bulk EDOS of the 
latter]^] In addition to just increasing the number of electrons absorbing the RF photons, such doping would 
also enhance the scattering intensity of excited electrons. Both of these tendencies are expected to result in 
enhanced HR. 




Energy relative to E p (eV) 


FIG. 13 Calculated spin-resolved local EDOSs against energy relative to Ep for Ta and Fe impurities in 


Au, reproduced from (Postnikov and Moldosanov 2012). Calculated EDOS of (Munoz et al. 2013) for 


Fe d-electrons in 32-atom quasirandom structure of Au 3 oFe 2 supports our results for Fe impurities in Au 
(Postnikov and Moldosanov 2012). 


3 Our calculations were done within the density functional theory, using the generalised gradient approximation for 
the exchange-correlation, by the SIESTA method (see http://departments.icmab.es/leem/siesta/ ), allowing full atom 
relaxation for the 64-atom supercell (4x4x4-replicated fee primitive cell with one atom substituted by impurity). 

4 The calculations done on 64-atom super cells (one impurity of Ta or Fe per 63 Au atoms) indicate that the total 
EDOS at the Fermi level is increased by ~50% compared to the pure Au. 
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Another imaginable way to bring about high EDOS would be due to /-electronic states residing at 
the Fermi level, that is the case of heavy-fermion compounds - see (Stewart 1984) for a review. Such 
compounds, typically possessing a rare-earth element as their essential constituent, exhibit an anomalously 
high residual resistivity at low temperatures, traceable to high electron scattering intensity, and, in its turn, 
- to a strong peak in the density of states of / electrons at the Fermi level. The presence of the compound 
CeCu6 among heavy-fermion systems permits to presume that gold, like a homologue of copper, may host 
/-element atoms so as to give rise to heavy-fermion behaviour. The nanoparticles of such compounds might 
possess a yet elevated heating rate. 


VI. DISCUSSION 


The analysis of conservation conditions for energy ( m e yAE e \ ss n vm - AE vm ) and momentum (i.e., for the 
movement direction of the excited electron after absorption of the RF photon and the LAVM) lets estimate 
the size of GNPs capable to produce high HR via a simultaneous involvement of several vibration modes 
with energies matching n vm -AE vnl , situated within the LAVMs of gold. It follows from our estimates that 
the optimal size of GNPs to use in RF hyperthermia is 5 — 7 nm. 

A number of known results seem to support our argumentation. 


(Moran et al. 2009), who explored 


heating of GNPs with sizes from 5 to 250 nm in the electric field of 13.56 MHz frequency, observed that the 
highest HR was inherent to GNPs of 5 nm size. ( Kruse et al.\ 2011) heated the GNPs by the 13.56 MHz 
electric field too, and, again, the smallest tested GNPs (~ 5 nm size) revealed the maximum HR. 

Other researches testify that as the GNPs sizes deviate from 5 — 7 nm, the HRs change from bad to 
worse. (Murdock and Hussain 2010) observed the heating of GNPs with the 4 nm and 13 nm diameters and 


pointed out that their HRs were several times lower than those for the 5 nm GNPs in the work by (Moran 
et al. |2009| ). (Li et al. 2011) exposed the GNPs with sizes of 20, 50 and 200 nm to the 13.56 MHz electric 
field and did not observe heating. In view of our analysis, the GNPs of 20 or 50 nm fall short of sizes (shown 
in Table [I]) that would enable a “direct” absorption of a RF photon; at the same time, these sizes are much 
larger than 5 — 7 nm, the “optimum” diameter for the phonon-assisted absorption. As for the 200 nm size, 
that falls close to a number appearing in Table |T| one can presume that the resonance for direct absorption 
is quite sharp and not satisfied in the real GNPs studied. 


(Liu et al. 2012) did not detect any significant heating of GNPs with sizes of 15 to 30 nm (i.e. larger 
than the optimum size) at the 13.56 MHz frequency. (Sirotkina et al. 2012) tried to treat tumours by the 
RF radiation of frequency 2.45 GHz using the gold nanorods of 30 nm diameter and 60 nm length, whereby 
no anticancer effect was identified. In our opinion, namely the large sizes of nanorods were the problem. 


(Corr et al. 2012a) and (Raoof et al. 2012) pointed out that heat generation was significantly reduced 


when the GNPs were aggregated. We attribute it to the fact that, as the GNPs come in contact, the common 
size of the formed aggregate is increased, and the condition ([5]) is not anymore fulfilled, resulting in reduction 
of HRs of GNPs. 

A valuable contribution to the analysis of experimental and theoretical studies in heating kinetics of 
GNPs in the RF range was offered by a recent work by (Letfullin et al. 2015). The authors have noted 


a spectacular failure of previous theoretical conclusions to account for the experimental evidence that RF 
radiation is capable of heating the GNPs. Our present work that emphasises the crucial role of LAVMs in 
the RF absorption suggests a plausible explanation of the above problem. 


(Huang et al. 2012) who studied in vivo the size dependence (within the 2 to 15 nm range) of the GNPs’ 
localization and penetration in cancer cells, multicellular spheroids and tumours, concluded that the 2 and 
6 nm GNPs demonstrated advantages over larger nanoparticles in terms of tumour uptake and permeability. 
Such GNPs were able to penetrate deeply into tumour tissue and provide high levels of accumulation in it. As 
an application of the smallest (~2 nm size) GNPs seems questionable due to loosening their biocompatibility 


(Leifert et al. 2013 Pan et al., 2012 2007), the 5 — 7 nm GNPs, possessing also the highest HRs, become 


very valuable instruments in the RF hyperthermia. 

Throughout the whole range of frequencies v used in medical or biological studies (^~10 MHz - 3 GHz), 
the RF photon energies hv are small compared to n vnl -AE vln of Eq. ([5]), for the LAVMs energy range of «3 
- 21 meV. Consequently, the condition © is generally satisfied, and the GNPs of about 5 — 7 nm size would 
have elevated HRs not only at 13.56 MHz, but everywhere within the mentioned interval of frequencies; in 
particular, at 2.45 GHz, the operation frequency of a domestic microwave oven with fine, inverter-regulated 
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power control. The use of the latter for laboratory investigations would permit performing genetic and 
biological experiments in vivo - e.g., with GNPs immersed in the tissues of larvae of such classical research 
target as Drosophila melanogaster. 


Recently, (San et al. 2013) reported that at the frequency of 0.35 MHz the Pt nanoparticles possessed 


HR twice higher than that of the GNPs and proposed to use them in RF hyperthermia. We attribute high 
HR of Pt nanoparticles to enhanced EDOS at the Fermi level of Pt, as compared to Au. For this reason, in 
Pt, a number of electrons able to absorb the RF photons is enhanced. 

One should expect that in conditions of the experimental capacitive electric transfer system used by 


(San et al. 2013) in nanoparticles of AuPd and AuPt alloys (Pd and Pt content within 40 - 70 mass %), 


the HRs would be even higher than those in Pt nanoparticles. This is because in these alloys, judging by 
their enhanced EDOS at Ep, the intensity of scattering of Fermi electrons is higher than that in pure Pt. 
However, an application of AuPd, AuPt, Pd and Pt nanoparticles in the RF hyperthermia is questionable 
due to their enhanced catalytic activities. 

Compared with these nanoparticles, the GNPs containing impurities of Ta and/or Fe, and hence devel¬ 
oping an enhanced EDOS at the Fermi level and, expectedly, an enhanced intensity of electron scattering, 


seem to be much more promising for the RF hyperthermia using the GNPs (Moldosanov and Postnikov 


2013 Postnikov and Moldosanov 2012). 


During last years, the use of GNPs as systems of drug and gene delivery into cancer cells has expanded 
considerably. An opportunity to synergistically combine these techniques with RF hyperthermia may open 
new terrains in contemporary cancer theranostics. Recently (Conde et al. 2015), apparently driven by 


a need to sense and overcome the cancer multidrug resistance, invented an implantable hydrogel with 
embedded DNA-coated GNPs of ~10 17 nm diameters. Our above discussion hints that an attempt to go 

towards slightly smaller particles, of 5 — 7 nm diameter, would “activate” the embedded GNPs for the RF 
hyperthermia, adding such an option to the GNPs’ other functions. 


Mirkin and co-workers (Chinen et al. 

2015 

Cutler et al. 

2012 

Jensen et al. 

2013 

Kouri et al. 2015[ 

Radovic-Moreno et al. 

2015 

Randeria et al. 

20151 developed and used spherical nucleic acid GNP con- 


Hypothetically, the 5 — 7 nm diameter GNPs along with optional RF heating thereof could be used in these 
cancer treating technologies, providing new research opportunities through temperature control. 

Another issue that enters the domain of feasible is the transfer of technologies developed for the plasmonic 


heating of GNPs onto the RF range. Zharov and co-workers (Galanzha et al. 2009 Nedosekin et al. 2013), 


aiming to prevent metastasis, used magnetic trapping of tumour cells circulating in the bloodstream with 
their simultaneous photoacoustic and photothermal detection. To this end, the gold-plated carbon nanotubes 
were employed. As it seems, the use of the 5 — 7 nm diameter GNPs instead of gold-plated carbon nanotubes 
may help to extend this method over applications in the RF range as well. 

A recent work by |Hemb ury et al.\ 2015) describes a promising so-called “quantum rattle”, that is, a 
hollow spherical particle (~150 nm of total diameter) with mesoporous silica shell (^25 nm thickness) host¬ 
ing both gold quantum dots (AuQDs) of < 2 nm diameter and GNPs (average crystallite size 7.3 nm in 
diameter). The quantum rattle is highly biocompatible and combines both cancer imaging and tumor treat¬ 
ment capabilities (chemotherapy and photothermal therapy in near infrared range). All these advantages 
in matching and even outperforming the state-of-the-art nanotechnology-based medical agents are achieved 
thanks to just AuQDs, whereas a contribution of GNPs is very small. Specifically, the photothermal therapy 
using the quantum rattles is provided through excitation of AuQDs by an infrared laser, then the AuQDs 
emit infrared fluorescence and enough heat to kill cancerous cells. In these quantum rattles, the GNPs are 
likely just processing waste. Meanwhile, as their sizes (7.3 nm) happen to fall close to the “favourable” 
range of 5 — 7 nm, the presence of GNPs in quantum rattles under discussion could be used for additional 
targeted heating of the latter by RF irradiation. 

Another recent work by (Tsai et al. 2015) outlines the use of gold nanorods (of 10 nm diameter / 41 nm 
length) connected to nanodiamonds, whereby the gold part serves to hyperthermia, by near-infrared laser 
heating, and the nitrogen-vacancy centers in diamond serve for temperature sensing. Here again, the use 
of nanoparticles of 5 - 7 nm size in place of nanorods would hopefully help to extend this technique over 
RF-induced hyperthermia. 

Last but not least, the RF heating of GNPs may be of interest for developers of therapeutic strategies 
targeting to inhibit amyloidogenic process in the Alzheimer’s disease. (Araya et al . 2008) suggested to use 
in this context the GNPs of 12.5±1.7 nm size heated by the microwave 14 GHz radiation. Chances are that 
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using the GNPs of “optimal” sizes would be useful for these tasks, too. 


VII. CONCLUSION 

We suggest a physical model of the size effect in heat generation in GNPs, which also accounts for 
a reduction of heat generation as the GNPs get aggregated. In this model, the LAVMs (dominating in 
the distribution of vibrational density of states) play an important role - an apparently novel element 
in the related theory framework. According to our model, the heating of GNPs is thought to consist of 
two consecutive processes: first, the Fermi electron absorbs simultaneously the RF photon and the LAVM 
available in the GNP; hereafter the excited electron is relaxed, exciting a LAVM with the energy higher than 
that of the previously absorbed LAVM. The model predicts that the GNPs to be effectively heated should 
possess diameters of ~ 5 — 7 nrn, i.e., very close to the experimentally inspected ~5 mn. The absorption 
band is expected to be very wide (~10 MHz - 3 GHz). This allows the use of frequencies typical for 
the “conventional” RF hyperthermia (without conducting nanoparticles). The energy release in the GNPs 
can be optimized by tuning the RF frequency, searching a compromise between the HR energy transfer 
efficiency and a penetration depth of the RF radiation into the biological tissue. The GNPs containing Ta 
or Fe impurity atoms are expected to be more effective heaters compared to nanoparticles of pure gold, due 
to enhanced electron density of states at the Fermi level. Gold nanoparticles with rare-earth impurity atoms 
are also brought into consideration as promising for the RF hyperthermia with conducting nanoparticles. 
The significance of the present study follows from the fact that the cancer specialists seeking approval for 
human clinical trials on the basis of their experimental results remain thus far in the dark in what regards 
the physical mechanism behind the observed trends. 
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Appendix A: Level spacing vs number of atoms in a GNP 


(Kubo 1962 1977) adopted the one-electron approximation to electrons in small metallic particles, 


whereupon D.Y. (Sun et al. 20011 predicted and J. (Sun et al. 2014) experimentally demonstrated the 


existence of the bulklike core region in them. This justifies considering the GNP as a crystalline one and 
permits to specify a relationship between the level spacing of electrons and number of atoms in GNP, using 
formulas for the free electron gas. A GNP of cuboid shape and consisting of N a atoms contains (one “free” 
electron per atom of gold) N a electrons filling the energy levels up to the Fermi energy E-p. The free-electron 
EDOS is (e.g., Wert and Thomson 1970): 


p(E) = dN'/dE = 47r(2m) 3/2 L 3 E 1/2 h 


-3 


(Al) 


where dN' is a number of states within the energy interval dE , L is the cuboid’s edge and h the Planck 
constant. Then the total number of electrons is 


Ep g 

N a = Jp(E)dE = y (2 mEpf/ 2 


(A2) 


In accordance with Eq. (Al), the spacing AE e \ between the Fermi level and nearest energy level above it 
equals ARj ss AN'/p{Ep) 1 where AN' = 2, hence 


all 


h 3 


Pi E f) 2n(2m) 3 / 2 L 3 E 3/2 ’ 


(A3) 
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From Eqs. (A2) and (A3), the Kubo’s formula is deduced: 

A-E e i 


4 Ep 
3 W a 


(A4) 


Appendix B: Counting vibration modes in GNPs of given size, within given energy range 

Assuming that the LAVMs propagate with (longitudinal) sound velocity ul along the circular contour 
at depth S under the surface of spherical particle of diameter E, the phonon energy quantum is AE vm = 
Vhh/[n(D — 2S)\. However, insofar as phonons are absorbed /created in the process of the electron excitations 
/ relaxations, the phonon energies can only change in blocks, commensurable with the electron states 
quantization, i.e., respecting the energy conservation condition of Eq. & ni e \^.E e \ ~ ?7. vm A£/ vm . Only 
certain combinations (ra e i,n V m) are possible, that selects the “resonant” values of D, as given by Eq. (ID- 
Table III lists the allowed groups (m e i, n vm , D) in the increasing order of D. 


Searching in Table III for the D values closest to the “reference” ones D i = 5 nm and D 2 = 10 nm, 
one finds D = 5.22 nm that comes along with n vm = 3, and D = 9.80 nm that comes along with n vm = 1. 
We note that the apparently “competitive” values of D by their closeness to the reference values, namely, 
D = 4.90 nm and 10.44 nm, can only be selected with much higher values of n vm and, consequently, may only 
intervene with much sparser distribution, and hence much smaller impact, or their allowed vibration modes. 
For D = 5.22 nm, n vm A£ vm = 3.35 meV, so there is no more than just one vibration mode, of the energy 
5x (n vm AE vm ), that falls within the FWHM range of the interest, 14.6 to 18.4 meV (see caption to Fig.[8]and 
the related text), hence = 1 in Eq. (19). For D = 9.80 nm, (n vm = 1)-A£ vm = 0.51 meV, and the 


above cited FWHM range may hosts much more modes at multipliers of this energy, namely, Ap^ HM f«7.5 
in Eq. ((191) . 
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